Formation of Liesegang patterns 
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It has been recently shown that precipitation bands characteristic of Liesegang patterns emerge 
from spinodal decomposition of reaction products in the wake of moving reaction fronts. This 
mechanism explains the geometric sequence of band positions a;„ ~ Q{1 + p)" and, furthermore, 
it yields a spacing coefficient p that is in agreement with the experimentally observed Matalon- 
Packter law. Here I examine the assumptions underlying this theory and discuss the choice of input 
parameters that leads to experimentally observable patterns. I also show that the so called width 
law relating the position and the width of the bands ui„ ~ Xn follows naturally from this theory. 
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I. INTRODUCTION 



Formation of precipitation patterns in the wake of 
moving reaction fronts (known as the Liesegang phe- 
nomenon) has been studied for more than a century |0-p| . 
The motivation for these studies has been diverse, coming 
from the importance of related practical problems such as 
crystal growth in gels, as well as from the fascination with 
a complex pattern that has eluded a clean-cut explana- 
tion (e.g. agate rocks are beheved to display Liesegang 
patterns). From a theoretical point of view, the main 
factor in the popularity was the belief that much can be 
learned about the details of precipitation processes (nu- 
cleation, growth, coagulation, etc.) by investigating the 
instabilities underlying this phenomenon. Currently, the 
Liesegang phenomenon is mainly studied as a nontrivial 
example of pattern formation in the wake of a moving 
front |Q-|| and there are speculations about the possibil- 
ity of creating complex mesoscopic structures using this 
rather inexpensive process. 

Liesegang patterns are easy to produce (Fig.l shows 
a particular experiment that we shall have in mind in 
the following discussion). The main ingredients are two 
chemicals A and B yielding a reaction product A -I- i? — > 
C that forms a nonsoluble precipitate C —> D under 
appropriate conditions [A — NaOH, B — MgCl2 and 
D = Mg{OH)2 in Fig.l]. The reagents are separated 
initially with one of them (_B, inner electrolyte) dissolved 
in a gel and placed in a test tube. Then at time t = Q 
an aqueous solution of the other reagent {A, outer elec- 
trolyte) is poured over the gel. The initial concentra- 
tion flo of A is chosen to be much larger than that of 
B (typically ao/6o ~ 10^), thus A diffuses into the gel 
and a reaction front moves down the tube. Behind the 
front, a series of stationary precipitation zones (Liesegang 
bands) appear at positions x„ (x„ is measured from 
the interface between the gel and the aqueous solution; 
n = 1,2, ...,10 — 20, typically). A band appears in a 
rather short time-interval thus the time of the appear- 
ance tn of the n-th band is also a well defined, experi- 



mentally measurable quantity. Finally, the widths of the 
bands Wn can also be determined in order to characterize 
the pattern in more detail. 
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FIG. 1. Liesegang patterns obtained with reagents 
A = NaOH and B = MgCh in polyvinylalcohol gel. The 
white precipitate is _D = Mg{OH)2- The height of the 
columns is 30cm and it takes about a 1-2 weeks for the pat- 
terns to form. The columns show different patterns due to 
the difference in the initial concentrations of the outer elec- 
trolyte NaOH. The experiments were carried out by M. 
Zrinyi (Technical University of Budapest). 

The experimentally measured quantities (a;„, t„ and 
Wn) in regular Liesegang patterns satisfy the following 
time-, spacing-, and width laws. 



Time law [Q: 



(1) 



This law is satisfied in all the experiments where it 
was measured and it appears to be a direct conse- 
quence of the diffusive dynamics of the reagents. 

• Spacing law [|| : The positions of the bands form a 
geometric series to a good approximation 



xn-Q{i+py 



(2) 



where p > is the spacing coefficient while Q is 
the amplitude of the spacing law. The quantita- 
tive experimental observations concern mainly this 
law. More detailed works go past the confirma- 
tion of the existence of the geometric series and 
study the dependence of the spacing coefficient on 



ao and bg. The results can be summarized in a rel- 
atively simple expression usually referred to as the 
Matalon-Packter law P,|lO[: 



p = F{bo) + G(5o) 






(3) 



where F and G are decreasing functions of their 
argument bo. 



Width law M 



Wn 



(4) 



This is the least established law since there are 
problems with both the definition and the mea- 
surement (fluctuations) of the width. Recent, good 
quality data ||l^ does support, however, the valid- 
ityof|). 

It should be clear that (|l|-|j) summarizes only those 
properties of Liesegang patterns that are common in a 
large number of experimental observations. There is a 
wealth of additional data on various details such as e.g. 
the secondary structures or the irregular band spacing 
[g|. These features, however, appear to be peculiarities 
of given systems. It is hard to characterize them and 
their reproducibility is often problematic as well. In view 
of this, it is not surprising that the theoretical expla- 
nations of Liesegang phenomena have been mainly con- 
cerned with the derivation of (13-0). 

The theoretical approaches to quasiperiodic precipi- 
tation have a long history and the two main lines of 
thoughts are called as pre- and post-nucleation theories 
(for a brief overview see |^). They all share the as- 
sumption that the precipitate appears as the system goes 
through some nucleation or coagulation thresholds. The 
differences are in the details of treating the intermediate 
steps "...C..." in the chain of reactions A + B -^ ...C... -^ 
D producing the precipitate D. In general, all the the- 
ories can explain the emergence of distinct bands but 
only the pre-nuclcation theories can account |p|,p|,[13|-[lq] 
for the time- and spacing laws of normal patterns. These 
theories are rather complicated, however, and have been 
developed only recently ||^ to a level that the dependence 
of p on the initial concentrations ao and &o can be investi- 
gated quantitatively, and connection can be made to the 
Matalon-Packter law [|,0. 

Unfortunately, there are several problems with the the- 
ories mentioned above. First, they employ a large num- 
ber of parameters and some of these parameters are hard 
to grasp theoretically and impossible to control experi- 
mentally (an example is the lower threshold in the density 
of C-s below which aggregation C + D ^ 2D ceases M). 
Second, some of the mechanisms invoked in the expla- 
nations are too detailed and tailored to a given system 



in contrast to the generality of the resulting pattern in 
diverse systems. A real drawback of the too detailed de- 
scription is that quantitative deductions are difficult to 
make even with the present computer power B . A final 
problem we should mention is the absence of an unam- 
biguous derivation of the width law in any of the theories. 
In order to avoid the above problems, we have recently 
developed a simple model of band formation Jl^l based 
on the assumption that the main ingredients of a macro- 
scopic description should be the presence of a moving 
reaction front and the phase separation that takes place 
behind the front. This theory contains a minimal number 
of parameters, it accounts for the spacing law, and it is 
simple enough that the existence of the Matalon-Packter 
law can be established numerically. The apparent success 
warrants a closer look at the model and, in this lecture, 
I will describe in detail how one arrives at such a model 
and what are the underlying assumptions of the theory. 
Then I would like to discuss the choice of input param- 
eters that yield experimentally observable patterns and, 
finally, I will show that the derivation of width law is 
straightforward in this theory. 



II. THE MODEL 



Let us begin building the model by taking a look at 
Fig. 1 . It shows alternating high- and low-density regions 
of the chemical Mg{OH)2 and the systems appear to be 
a quasi-steady state (actually, there are experiments that 
suggest that the pattern does not change over a 30 years 
period H ) . We shall take this picture as an evidence that 
phase separation |I^ underlies the formation of bands 
and, furthermore, that the phase separation takes place 
at a very low effective temperature (no coarsening is ob- 
served). 

The phase separation, of course, must be preceded by 
the production of C-s. This is the least understood part 
of the process and it is particular to each system. What 
is clear is that due to the condition oq S> &o a reaction 
front {A -\- B ^> something) moves down the tube diffu- 
sively (note that this is the point where the role of the gel 
is important since it prevents convective motion). The 
result of the reaction may be rather complex (intermedi- 
ate products, sol formation, etc.) and one of our main 
assumptions is that all these are irrelevant details on a 
macroscopic level. Accordingly, the production of C will 
be assumed to be describable by the simplest reaction 
scheme A + B ^ C. 

Once A + B ^' C is assumed, the properties of the 
front and the production of C-s are known Jl^ . Namely, 
the front moves diffusively with its position given by 
Xf = y/2Dft, the production of C-s is restricted to a 
slowly widening narrow interval [wf{t) — wot^'^] around 



Xf, and the rate of production S{x,t) of C-s can be ap- 
proximated by a gaussian (the actual form is not a gaus- 
sian, see pa] for details about a non-moving front) 



S{x,t) 



So 



exp 



^f{t)? 



2wyt) 



(5) 



The parameter Df can be expressed through ao, bo, and 
the diffusion coefhcients of the reagents {Da,Di,) while 
So and wo depends also on the rate constant, k, of the 
reaction A + B ^ C. 

An important property of the front is that it leaves 
behind a constant density cq of C-s |g] and cq depends 
only on oq, bo, Da and Di,. This is important because the 
relevant parameters in the phase separation are Df and 
Co (where and how much of the C-s are produced |20|) 
and thus the least available parameter (fc) does not play 
a significant role in the pattern formation. 

Having a description of the production of C-s, we must 
now turn to the dynamics of their phase separation. Since 
the emerging pattern is macroscopic, we shall assume 
that, on a coarse-grained level, the phase separation can 
be described by the simplest 'hydrodynamical' equation 
that respects the conservation of C-s. This is the Cahn- 
Hilliard equation |2^ or, in other context, it is the equa- 
tion for model B in critical dynamics |22| . This equation, 
however, requires the knowledge of the free-energy den- 
sity {T) of the system. For a homogeneous system, JF 
must have two minima corresponding to the low- (c; ) and 
high-density (c/j) states being in equilibrium (Fig.l). The 
simplest form of J^ having this property and containing a 
minimal number of parameters is the Landau-Ginzburg 
free energy (Fig. 2) 



J- = —-em?' + —jrn + -cr(Vm)^ , 



(6) 



where in = c — {ci + Ch)/2 is the density, c, of the C-s 
measured from the average of the two steady state values 
(we are following the notation in pq] where the 'mag- 
netic language' has its origin in a connection to Ising 
lattice gases). The parameters e, 7, and a are sys- 
tem dependent with e > ensuring that the system is 
in the phase-separating regime, a > provides stabil- 
ity against short-wavelength fluctuations, and requiring 
\/e/'j = (c/i — q)/2 fixes the minima of J- at ztnie corre- 
sponding to Q and Ch- Note that the m —^ —m symmetry 
is usually not present in a real system and J-" could con- 
tain e.g. an rrr' term. The presence or absence of the 
m —^ —m symmetry, however, is not relevant for the dis- 
cussion that follows. 
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FIG. 2. The homogeneous part of the free energy as a func- 
tion oi m = c— {ci + Ch)/2. The phase separation is an acti- 
vated process in the metastable regimes while it goes by spin- 
odal decomposition in the (linearly) unstable regime. The 
spinodal point (o) separates these regimes. The arrows are 
meant to illustrate how the density in the front increases to- 
wards mo and how the phase separation to the steady states 
(zhrrie) takes place when the density reaches the spinodal value 



Using (g) and including the source term, the Cahn- 
Hilliard equation takes the form 



dtm 



-AA ( em — 7TO -I- a Am ) + S* . 



(7) 



were A is a kinetic coefficient. The above equation should 
contain two noise terms. One of them should be the ther- 
mal noise while the other should originate in the chemi- 
cal reaction that creates the source term. Both of these 
noise terms are omitted here. The reason for neglecting 
the thermal noise is the low effective temperature of the 
phase separation as discussed in connection with Fig.l. 
The noise in S, on the other hand, is dropped since the 
A + B ^ C type reaction fronts have been shown to be 
mean- field like above dimension two |23|| . 

The absence of noise means that the phase separation 
can occur only through spinodal decomposition |0,e3. 
Thus the assumption behind omitting the noises is that 
the characteristic time of nucleation is much larger than 
the time needed by the front to increase the density of 
C-s beyond the spinodal value {—mg in Fig. 2) where the 
system is unstable against linear perturbations. Since 
there are examples where the bands appear to be formed 
by nucleation and growth H, the spinodal decomposi- 
tion scenario is clearly not universally applicable, and 
one should explore the effects of including noise (this be- 
comes, however, an order of magnitude harder problem). 

Eq. (R) together with the form of the source (H) defines 
now our model p6| that produces regular Liesegang pat- 
terns (Fig. 3) satisfying the spacing law (0) and, further- 
more, the spacing coefficient is in agreement with the 
Matalon-Packter law (ph. Fig. 3 shows a rather general 



picture that is instructive in understanding the pattern 
formation. The last band acts as a sink for neighbor- 
ing particles above —me (cj) density. Thus the C-s pro- 
duced in the front end up increasing the width of the last 
band. This continues until the front moves far enough so 
that the density in it reaches the spinodal value. Then 
the spinodal instability sets in and a new band appears. 
Remarkably, the above picture is rather similar to the 
phenomenological 'nucleation and growth' scenario Iq] 
with the density at the spinodal point playing the role 
of threshold density for nucleation. It is thus not entirely 
surprising that both of these theories do equally well in 
producing the spacing- and the Matalon-Packter law. 

One should note that the actual form of ^ does not 
play an important role in the picture developed above. 
The crucial feature is the existence of a spinodal density 
above which phase separation occurs. This is the mean- 
ing of our previous remark about the irrelevance of the 
rrt' term in the free energy (of course, one should also 
realize that explanations of details in experiments may 
require the inclusion of such terms). 
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FIG. 3. Liesegang pattern obtained for front parameters 
Df = 21.72, Wo = 4.54, and ^o = 0.181 with length, time, 
and m (concentration) measured in units of \/ o je — 2-10~''m, 

a/(Ae^) — 40s, and ^J e/'y, respectively. The dashed line de- 
notes the rate of production of C-s (S), measured in units of 
\e^^'^/(')^''^a) and magnified by a factor 2 ■ lO''. The dotted 
line represents the density at spinodal point, —nis = — 1/\/3. 



III. CHOICE OF PARAMETERS 

Fig. 3 shows the results of numerical solution of eq.ffl) 
with the same parameter values as in Fig. 3 of ref. pq ] 
but stopped at an earlier time so that the visual similar- 
ity to the experiments (number of bands in Fig.l) would 
be greater. In this section, we shall examine whether the 
parameters used for obtaining this resemblance have any 
relevance to real Liesegang phenomena. 



The experimental patterns have a total length of about 
^exp ~ 0.2m and the time of producing such a pattern is 
about 1-2 weeks (we shall take Texp ~ lO^s). Since our 
model has a length-scale Ith = ^/cr/e and a time-scale 
Tth — o'/(Ae^), they can be chosen so \\J o je = 2 • 10~^m 
and cr/(Ae^) = 40s] that i^xp ~ tth and r^xp ~ Tth- Once 
we have chosen tth and Tth we can start to calculate other 
quantities and see if they have reasonable values. 

It is clear from Fig. 3 that the widths of the bands are 
in agreement with the experiments, they are of the order 
of a few mm at the beginning and approach to ^ 1cm at 
the end. The width of the front is also of the order of 
1cm after lO^s. Unfortunately, there is no information on 
the reaction zone in this system. In a study of a differ- 
ent system |25| ] it was found that ■u;/(t=2 hours) w 2mm. 
Extrapolating this result tot = lO^s one finds Wf « 1cm 
[note that the exponent of the increase of Wf{t) ^ t^/^ is 
small] in agreement with the observed value. 

Next we calculate the diffusion coefficient of the front, 
Df = 21.72 • fth/^^th « 2 • 10-^mVs. This value ap- 



'f 



pears to be an order of magnitude larger than the usual 
ionic diffusion coefficients (D « 10~^m^/s). One should 
remember, however, that Fig. 3 is the result for initial 
conditions ao/bo = 10^ [Q and, for this ratio of ao/bo, 
the diffusion coefficient of the front Df is about 10 times 
larger than Da {Df/Da « 10 see Fig.4 in §). Thus Df 
also comes out to be the right order of magnitude. 

We do not have information on the amplitude (5*0) 
of the source but, once the concentrations (aq, 60) are 
given and Df and Wf are known then Sq is fixed by the 
conservation law for the C-s. Thus the correct order of 
magnitude for Df and Wf should ensure that 6*0 is also 
of right order of magnitude. 
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FIG. 4. Details of the time-evolution of the formation of 
the last band in Fig. 3. The time to is the moment when the 
concentration at the front reached the spinodal value and St 
is measured from to- 

Finally, we shall calculate the time it takes for a band 



to form. It is well known that the bands appear rather 
quickly. From the visual notice of the beginning of the 
band formation it takes about Ti„i — 30 — 60 minutes 
for the band to be clearly seen and then it takes much 
longer to increase its width to the final value. In order 



{Xn+l - Xn)co = [Xn+l - Xn - Wn)ci + WnCh 



(8) 



to calculate 



let us consider the formation of the last 



band in Fig. 3 (see Fig. 4). The lower limit of the den- 
sity that can be visually noticed is, of course, not well 
defined. We shall assume that this density corresponds 
to TO = i.e. it is the halfway density from q to c^. 
This means that we see the beginnings of the band at 
5t = 410min and the density reaches well above 90% of 
its final value by 6t = 470min. Consequently, we ob- 
tain again an estimate {jini ~ 60inin) for an observed 
quantity that agrees with the experiments. As a result of 
the above estimates, we feel that the parameters in our 
model can indeed be chosen so that they are relevant to 
real Liesegang experiments. 



IV. WIDTH LAW 

The width law is problematic from experimental point 
of view since the fluctuations in the widths appear to be 
large. Part of the difficulties are undoubtedly due to the 
fact that the boundaries of the bands are not sharply de- 
fined and high-resolution digitizing methods are needed 
in a precise analysis. The most thorough experiment to 
date has been carried out recently lia] with the result 
Wn ~ a:" where a ~ 0.9 — 1.0. 

As to the theories, they also have their share of difficul- 
ties since, on a microscopic level, the growth of the width 
involves precipitation processes in the presence of large 
concentration gradients, while a macroscopic treatment 
must elaborate on the dynamics of the interfaces between 
two phases. Accordingly, there are only a few works to 
report on. Dee IJ] used reaction-diffusion equations sup- 
plemented by terms coming from nucleation and growth 
processes and obtained w„ ^ x„ from a rather limited 
(6 bands) numerical result. Chopard et al. g employed 
cellular automata simulations of a phenomenological ver- 
sion of the microscopic processes and found w„ ~ a:" with 
a ~ 0.5 — 0.6. Finally, Droz et al. [Q combined scaling 
considerations with the conservation law for the number 
of C particles to obtain a in terms of the scaling proper- 
ties of the density of precipitates in the bands. Assuming 
constant density they found a = 1. Our derivation below 
parallels this last work in that the same conservation law 
is one of the main ingredient in it. 

In our theory, the derivation of the width law is 
straightforward. One combines the facts that (i) the re- 
action front leaves behind a constant density (cq) of C-s, 
(ii) the C-s segregate into low (q) and high (c/i) density 
bands, (iii) the number of C-s is conserved in the segre- 
gation process; and writes down the equation expressing 
the conservation of C-s 



Using now the spacing law (g) that has been established 
for this model one finds 



_ p(co - q) 
Ch - ci 



CyXr, 



(9) 



We have thus derived the width law and obtained the 
coefficient of proportionality, C, as well. The importance 
of C lies in that measuring it provides a way of accessing 
cg that is not easily measured otherwise. 



V. FINAL REMARKS 

In summary, we have seen that the spinodal decom- 
position scenario for the formation of Liesegang patterns 
performs well whenever quantitative comparison with ex- 
periments is possible. It remains to be seen if the applica- 
bility of this model extends beyond the re(7w/ar patterns. 
One should certainly try to use this theory to explain 
the exotic patterns (e.g. inverse patterns, helixes) that 
are experimentally reproducible and lack even qualitative 
understanding. 
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